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Abstract. The S^-Tools are a set of Python-based routines and interfaces whose purpose is to provide user- 
friendly access to the SKA Simulated Skies (S ) set of simulations, an effort led by the University of Oxford in the 
framework of the European Union's SKADS program (http://www.skads-eu.org). The databases built from the 
S simulations are hosted by the Oxford e- Research Center (OeRC), and can be accessed through a web portal at 
http://s-cubed.physics.ox.ac.uk. This paper focuses on the practical steps involved to make radio images from 
the 5 3 -SEX and 5 3 -SAX simulations using the 5' 3 -Map tool and should be taken as a broad overview. For a more 
complete description, the interested reader should look up the user's guide. The output images can then be used as 
input to instrument simulators, e.g. to assess technical designs and observational strategies for the SKA and SKA 
pathfinders. 
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1. Introduction 

The SKA Simulated Skies project (S 3 ) aims at building 
mock radio skies suitable for planning science with the 
SKA and the many pathfinder experiments. The project 
comprises a number of simulations, of which essentially 
two are relevant for this paper: 

• 5 3 -SEX (Semi-Empirical cXtragalactic) : a large- 
scale simulation of the extragalactic radio continuum sky 
covering a sky area of 20° x 20°, out to redshift z = 20, 
and down to 10 nJy It also includes Hi gas masses 
for star-forming galaxies. For a complete description, see 



Wilman et al. (2008) 



• S 3 -SAX (Semi- Analytical eXtragalactic) : a smaller- 
scale simulation of the extragalactic Hi and CO 
line emissions, derived from the Millenium simulation 
( |Springel et al. (2005) [ ), which comes in two flavours. S - 
SAX-Sky is a skyfield simulation giving the apparent 
properties of galaxies, with a field of view that depends 
on the maximum redshift, while S^-SAX-Box is a sim- 
ulation of a cubic volume giving their intrinsic proper- 
ties. In all of this paper, 5 3 -SAX shall be understood 
as meaning S^-SAX-Sky. For a complete description, see 



Obreschkow et al. (2009) 



2. Querying the S 3 databases 

2.1. Database structure 

The databases can be queried via SQL forms on the 
dedicated S 3 -SEX and S 3 -SAX sections of the S 3 web- 
site hosted at |http: / / s-c ubed .physic s . ox. ac .uk} The 
queries should be written in the SQL format, and query 
examples are given. 

The S 3 -SEX and S 3 -SAX databases have similar struc- 
tures, and consist of respectively three and two tables. 



• For 5 3 -SEX : a Clusters table holding cluster prop- 
erties, a Galaxies table holding galaxy properties (possibly 
including the cluster index to which they belong) and a 
Components table, which lists the properties of the sev- 
eral components that may make up a single galaxy, such 
as cores and radio lobes. 

• For 5 3 -SAX : a galaxiesjine table containing the ap- 
parent position and emission line properties of the galaxies 
and a galaxies_delu table containing the intrinsic properties 



of the de Lucia et al. (2006) catalog 



For more information on the structure of the database 
and the attributes held in the various tables, please refer to 
the appropriate sections of the S 3 website. It is advisable 
to do so before reading on, as in the following we shall 
mention some of these attributes explicitly. 

2.2. Mandatory attributes 

To make maps from these simulations, it is mandatory 
that some properties - listed below - be included in the 
query, as the mapping algorithms require them : 

• For S®-SEX : 

In the Galaxies table : galaxy, redshift or modified_redshift, 
distance, sftype, agntype (these latter two are necessary if 
one wishes to make separate images for the different types 
of sources), m_hi (if Hi mapping is to be performed). 
In the Components table : galaxy, right_ascension, decli- 
nation, position_angle, major_axis, minor_axis, L151, L610, 
L1400, L4860, L14000. 

• For S 3 -SAX : 

In the Galaxies table only : ra, deel, zapparent, diskpo- 
sitionangle, diskinclination, rmolc, distance. To these, one 
should add hiintflux, himajoraxis_10max, hiwidth50 and hi- 
width20 for mapping the Hi line, as well as cointflux_J, 
h2majoraxis_10max, cowidth50 and cowidth20 for mapping 
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the C0( J — * J — 1) line. When using the Kapteyn Hi tem- 
plates (see below), hubbletype should also be retrieved. 

2.3. Structure of query result files 

Query results are saved in a gzipped tarball id.tar.gz, 
where id is a hash given by the server at submission 
time. The tarball contains two plain text files : id. sql, in 
which the original SQL-formatted query is recalled, and 
id. result, which contains the comma-separated query 
results proper. 

3. Installing the S* 3 -Tools 

3.1. Prerequisites and installation 

The £ 3 -Tools are currently hosted on an external website 



S3Map 



at http : / / www . Ira . ens . f r/ ~levrier/Recherche/S3/ 



although they should find a permanent home on the S' 
webserver in the near future. 

To install the 5 3 -Tools, you should refer to that web- 
site or to the README file included in the distribution. 
In particular, the user should be aware that the S^-Tools 
require a number of python libraries, namely the scipy, 
numpy, os, math and pyfits packages. For Mac users, these 
are all included in the SciSoft package (at least from ver- 
sion 2008.3.1 onwards) for OS X. 

Provided these packages are functioning properly on 
your system, installation is a simple matter of untar- 
ring the S3Tools . tgz archive and setting a few paths 
in the Conf ig_Path.py file in the Config/ subdirec- 
tory. In the following, we shall assume that instal- 
lation is done in a directory named S3Tools under 
your home directory. Configuration files should there- 
fore be in ^VS3Tools/Conf ig/ and the main routines in 
~/S3Tools/Routines/. 

3.2. Templates 

To make line maps from S^-SAX, Hi and CO templates 
are required. These were made by D. Obreschkow and are 
available on the same webpage. Be warned that the archive 
weighs over 2 GB. A different set of templates for Hi has 
been produced by R. Boomsma from Kapteyn Institute 
and are also readily available to use. 

3.3. Global Sky Model 

The Global Sky Model (GSM) - also dubbed S 3 - 
GAL on the S* 3 website - is a model of the 
Galactic foregrounds from 10 MHz to 100 GHz by 
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Fig. 1. The S" 3 -Map graphical user interface 

installation directory specified via the GSMJDIR variable 
in Conf ig_Path.py. 

4. Making maps and spectral data cubes 

4.1. The S 3 -Map GUI and underlying scripts 

Making maps from query results is done via S^-Map, a 
graphical user interface (GUI) written in python, which is 
run by moving to the Routines/ subdirectory and typing 
python S3Map.py at the command line. A window simi- 
lar to the one presented in Fig. [1] should appear after a 
few moments. It is through this interface that users spec- 
ify which maps to make. To be more specific, the state of 
the GUI is passed as arguments to a python script com- 
mand, either SEX_map_script . py or SAX_map_script . py. 
It is this command that is in charge of actually building 
the maps using lower-level routines, and it may be printed 
out on the screen by clicking on the Show button. The user 
can then copy-paste this command into a terminal window 
and perform the same mapping non-interactively. This is 
a useful trick to know for building batch jobs or debugging 
purposes. For a lengthy description of the syntax used for 
the arguments, please refer to the S^-Tools user's guide. 



deOliveira-Costaet al. (2008)| It is possible to in- 4 2 . Common workings of the S 3 -Map GUI 



elude these foregrounds in spectral cubes built from 
the SKA Simulated Skies with S 3 -Map. To be 
able to use this option, the GSM data and rou- 
tines should be installed from the GSM website 
http://space.mit.edu/home/angelica/gsm/, and the 



Click on the Browse button to select a tarball containing 
results from a database query. This prompts «S 3 -Map to 
check which database was queried, so that some of the 
buttons, menus, entries and options become either active 
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or inactive, depending on their relevance. 5 3 -Map also fills 
in the mapping area fields with extremal values for right 
ascension (X) , declination (Y) and redshift (z) found in the 
query results. The From Query checkboxes are then set, 
but the user may deactivate them and manually specify a 
different mapping area. 

Regarding frequency configuration, this is done via 
three entry fields : Minimum Frequency (central frequency 
of the first channel), Frequency Resolution (full channel 
width) and the number of channels. For line mapping, a 
convenient Setup button has been included : Simply se- 
lect which line in which redshift range to look for and the 
number of channels to use, Setup automatically adjusts 
Minimum Frequency and Frequency Resolution for you. 

The user also has the possibility of putting in the ef- 
fects of a Gaussian beam, via the Beam Major FWHM, 
Beam Minor FWHM and Beam Position Angle entry fields. 
Please note that a beam smaller than the pixel size spec- 
ified in the corresponding field has of course no effect. In 
that case, no convolution is performed. 

Gaussian noise can be added to each frequency plane 
in the output spectral cube. This is done via the dedicated 
Noise level entry field. Currently, the said noise level sets 
the ratio of the added noise's rms value to the maximal 
value of the noiseless map. 

It is also possible to add the Galactic foreground signal 
at each frequency form the Global Sky Model described in 
13.31 For this, simply tick the Global Sky Model checkbox. 
Be aware that this is much lower spatial resolution than 
the SKA Simulated Skies, and so may appear in most cases 
as just an added constant over the requested skyfield. 

The unit of the third axis can be set, using the Z- 
axis menu, to one of frequency, velocity or wavelength. 
The user should be aware that the "velocity" option only 
applies to the mapping of a single line. If this option is 
selected while requesting the mapping of several lines or 
no line, the script will revert to a frequency axis. 

The usual convention for the orientation of the right 
ascension axis is that it increases to the left (i.e. to the 
east). This can be changed by ticking the Reverse R.A. 
checkbox. 

Finally, to use the mapped field as input to instru- 
ment simulators, the user has a possibility to specify the 
pointing corresponding to the center of the map. This is 
done using the Map Center entry fields at the bottom of 
the GUI. When the checkbox next to it is off, this point- 
ing is not used and the positions found in the query result 
are used. Considering that the simulations are centered on 
(0,0), this may lead to strange behaviour in later stages 
of the end-to-end simulations, and it is therefore not rec- 
ommended. 



Fig. 2. Dealing with image accuracy for few-pixel compo- 
nents. Left : without regridding ; Right : with regridding. 

4.3. Some points specific to S 3 -SEX 
4.3.1. Elliptical templates 

In S^-SEX, sources are made up of components that are ei- 
ther ellipses or point sources, and each line in the query re- 
sults file corresponds to such a component. As mentioned 
in !2.2[ the information given for each component includes 
position, major and minor axes, position angle, continuum 
fluxes at the reference frequencies, and possibly Hi mass. 
Given a window and a pixel size for the desired map, we 
can build a template image for each single component by 
setting to one the values of pixels whose centers fall within 
the ellipse, and to zero all other pixel values. 

This approach is fine if the component's size is much 
smaller or much larger than a pixel. In the former case only 
one pixel is non-zero, and in the latter case, the compo- 
nent's border in the template image resembles the desired 
ellipse quite closely. A problem arises with components 
whose sizes are a few pixels across, as can be seen on the 
left side of Fig. H 

To correct for this, the grid is temporarily refined, as 
on the right side of Fig. [21 with a minimum ratio of semi- 
major axis to refined pixel size given by the value of the 
Image accuracy entry field on the GUI. Then each coarse 
pixel in the original grid is given a value proportional to 
the total number of non-zero refined pixels that belong to 
this coarse pixel. 

In all cases, the template image is then normalized so 
that the integrated flux over the component is 1 Jy, and 
convolved by a flux-conserving beam. 



4.3.2. Continuum and Hi fluxes 

Simultaneously, the continuum flux at the desired fre- 
quency - or frequencies - is computed from a log- 
polynomial interpolation based on the reference fluxes re- 
trieved from the query results. As for the total Hi flux, it 
has to be derived from the Hi gas mass and distance using 
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"That Setup button is a stroke of genius" (I. Hey wood 
dixit) . 



and a synthetic double-peaked profile, is built randomly 
using some prescriptions based on galaxy typH The Hi 
flux in each frequency channel is computed by integration 



For more detail, please refer to the user's guide. 
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of this synthetic profile over the channel's width. The to- 
tal (continuum+Hi) flux is then used to multiply the tem- 
plate image at each of the desired frequencies, to obtain a 
spectral data cube which is finally pasted onto the overall 
cube. This means that the emission profile is painted uni- 
formly over the galaxy's disc. This approach is sufficient 
at low spatial resolution. 

4.3.3. Polarized continuum emission 

By default, only total intensity continuum emission is 
mapped, but in the l,Q,U,V menu, the user may specify 
that linearly polarized emission should be mapped as well. 
This emission is computed from the source's total intensity 
and observing frequencies by routines included in the S 3 - 
Tools and which were kindly provided by Jorn Geisbiisch 
(Cavendish Laboratory, Cambridge University). In that 
case, the output cube has a fourth dimension, with sepa- 
rate Stokes I, Q and U hyperplanes. In the algorithm, for 
each source, the same template built for total intensity is 
used for Stokes Q and U, it is simply multiplied by the 
appropriate fluxes. 

4.3.4. Source types 

There are 5 different source types in the S^-SEX simula- 
tion (Radio-quiet AGN, FRI and FRII radio-loud sources, 
quiescent star-forming galaxies and starbursting galax- 
ies) which are distinguished via their sftype and agn- 
type attributes. The user can choose to map only se- 
lected types by using the Source Types pull-down menu. 
Obviously, only those source types which were queried can 
be mapped, and the sftype and agntype must have been 
included in the query for the selective mapping to be pos- 
sible. On a related note, GPS sources are a subclass of 



the FRI and FRII source types (see Wilman et al. (2008) ) 



The user may specify whether these should be included or 
excluded, using the GPS menu. 

4.4. Some points specific to S 3 -SAX 

The Hi and CO emission maps from »S 3 -SAX data require, 
as mentioned in !3.2l a number of position-position- velocity 
(PPV) templates. For Hi, the user has the choice between 
templates from D. Obreschkow (Oxford option in the HI 
Data menu) or templates from R. Boomsma (Kapteyn op- 
tion). For CO, only templates by D. Obreschkow are avail- 
able. They are included in the archive along with the Hi 
templates. 

Ten J — > J — 1 rotational transitions of CO (J go- 
ing from 1 to 10) are available in the database and can 
be queried. The queried lines are active in the Transitions 
menu, and by default they are selected for mapping. 

In all cases, for a given galaxy, the relevant template is 
retrieved based on the values of hubbletype, rmolc, himajo- 
raxis_10max, h2majoraxis_10max, diskinclination, cowidth50 
and hiwidth50, then submitted to necessary up- or down- 



scaling in the spatial and velocity dimensions, rotation and 
normalization to the integrated line flux. The template is 
eventually pasted onto the output cube. 

5. Upcoming extensions 

The remaining items on the agenda for the S^-Tools are : 

• Inclusion of the epoch of reionization (EoR) signals 
produced independently by M. Santos (1ST Lisbon) and 
B. Semelin (Observatoire de Paris / LERMA). This will 
be very similar to the way the GSM is included. 

• Inclusion of polarized emission in S* 3 -SAX via tem- 
plates produced by R. Beck and T. Arshakian (MPIfR) 

• The mapping scripts currently deal with one source 
at a time, building or finding a template, transforming it 
and pasting it onto the master cube. Obviously, this could 
be parallelized. 

• An older version of the S 3 -Tools included annotation 
files for the KARMA visualisation software. This should 
eventually be re-included. 

• The mapping scripts should be installed on the S 3 
web server and a front-end form similar to the GUI de- 
veloped for the website. This would allow users to request 
maps from the server without having to download and 
install the S 3 -Tools. 

• A better treatment of Hi in 5 3 -SEX is possible, for 
instance by using Hi surveys such as HIPASS to derive 
synthetic line profiles consistent with observations. 

Acknowledgements. This work was supported by the European 
Commission Framework Program 6, Project SKADS, Square 
Kilometre Array Design Studies (SKADS), contract no 011938. 

References 

de Lucia, G., Springel, V., White, S.D.M., Croton, D., 
Kauffmann, G., 2006, MNRAS, 366, 499 

de Oliveira-Costa, A., Tegmark, M., Gaensler, B.M., Jonas, J., 
Landecker, T.L., Reich, P., 2008, MNRAS, 388, 247 

Obreschkow, D., Klockner, H.-R., Heywood, I., Levrier, F., 
Rawlings, S., 2009, ApJ, 703, 1890 

Springel, V., White, S.D.M., Jenkins, A., Frenk, C.S., Yoshida, 
N., Gao, L., Navarro, J., Thacker, R., Croton, D., Helly, 
J., Peacock, J. A., Cole, S., Thomas, P., Couchman, H., 
Evrard, A., Colberg, J., Pearce, F., 2005, Nature, 435, 629 

Wilman, R.J., Miller, L., Jarvis, M.L., Mauch, T., Levrier, F., 
Abdalla, F.B., Rawlings, S., Klockner, H.-R., Obreschkow, 
D., Olteanu, D., Young, S., 2008, MNRAS, 388, 1335 



